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Abstract 



Singularly perturbed vector nonlinear Schrodinger equations (PVNLS) are inves- 
tigated. Emphasis is placed upon the relation with their restriction: The singularly 
perturbed scalar nonlinear Schrodinger equation (PNLS) studied in [10]. It turns 
out that the persistent homoclinic orbit for the PNLS [10] is the only one for the 
PVNLS, asymptotic to the same saddle. 
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1 Introduction 



In recent years, novel results have been obtained on the solutions of the vector nonlinear 
Schrodinger equations [3] [4] [17]. Abundant ordinary integrable results have been carried 
through [16] [6] , including linear stability calculations [5] . Specifically, the vector nonlinear 
Schrodinger equations can be written as 

ipt + Pxx + ^Ol 2 + x\q\ 2 )p = o, 

iqt + qxx + ^{x\p\ 2 + \q\ 2 )q = 0, 

where p and q are complex valued functions of the two real variables t and x, and \ 
is a positive constant. These equations describe the evolution of two orthogonal pulse 
envelopes in birefringent optical fibers [14] [15], with industrial applications in fiber com- 
munication systems [7] and all-optical switching devices [8]. For linearly birefringent fibers 
[14], x — 2/3. For elliptically birefringent fibers, x can take other positive values [15]. 
When x = 1, these equations are first shown to be integrable by S. Manakov [13], and 
thus called Manakov equations. When x is not 1 or 0, these equations are non-integrable. 
Propelled by the industrial applications, extensive mathematical studies on the vector 
nonlinear Schrodinger equations have been carried. This article is another piece of math- 
ematical works. Like the scalar nonlinear Schrodinger equation, the vector nonlinear 
Schrodinger equations also possess figure eight structures in their phase space. Figure 
eight structures are also called separatrices. For two-dimensional Hamiltonian systems, 
figure eight structures are given by the singular level sets of the Hamiltonian. Our goal 
is to understand the consequence of such figure eight structures when these equations are 
under perturbations, in particular, whether or not the consequence is chaotic dynamics. 
Specifically, we will investigate the type of chaos created through homoclinic orbits, (cf: 
the section "Conclusion and Discussion".) 

This article will deal with the problem on the existence of homoclinic orbits for the 
singularly perturbed vector nonlinear Schrodinger equations, 

ipt + Pxx + \{\p\ 2 + \q\ 2 )p = ie[pxx - ap - pe^} , (1.1) 

iqt + qxx + \{\p? + \q\ 2 )q = ie[q xx - aq - , (1.2) 

where p(t, x) and q(t, x) are subject to periodic boundary condition of period 2%, and are 
even in x, i.e. 

p(t, x + 2tt) = p(t, x) , p(t, -x) = p(t, x) , 

q(t, x + 2tt) = q(t, x) , q(t, -x) = q(t, x) , 

uj G (1, 2), a > and j3 are real constants, and e > is the perturbation parameter. 
By the transformation, 

p = pe l ^ 2t , q = qe^ 2t , 
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and dropping , one gets an autonomous system, 

ipt+Pxx + ^[(bl 2 + |<?| 2 ) - u 2 ]p = ie[p xx -ap-fi] , (1.3) 

iqt + q xx + ^[{\p\ 2 + \q\ 2 ) - ^ 2 }q = iefoxx -aq- (5\ , (1.4) 

The rest of this article will be dealing with this form of perturbed vector nonlinear 
Schrodinger equations (PVNLS). 

The corresponding work for regularly perturbed scalar nonlinear Schrodinger equation 
has been accomplished in [11], where the singular perturbation ed 2 is replaced by a regular 
perturbation ed 2 which is a bounded Fourier multiplier obtained from a truncation of 
ed 2 . The main difficulty introduced by the singular perturbation ed 2 x is that it breaks the 
spectral gap condition of the unperturbed system. Thus the standard invariant manifold 
results do not hold. Nevertheless, certain invariant manifold results can be established, 
which are sufficient for establishing the existence of homoclinic orbits. The corresponding 
work for such singularly perturbed scalar nonlinear Schrodinger equation (PNLS) has 
been accomplished in [10]. The crucial new features of this article are listed below: 

1. Homoclinic orbits asymptotic to the spatially uniform and temporally periodic 
solutions up to phase shifts, are constructed for the integrable vector nonlinear 
Schrodinger equations ((1.1, 1.2) at e = 0). In particular, such homoclinic orbits 
are built with spatially non-periodic and non-anti-periodic Bloch eigenfunctions. 
This study provides an elegant formula which is missing in [16] [6] [5]. 

2. For the scalar nonlinear Schrodinger equation, a complete sequence of Melnikov 
vectors can be constructed out of the Floquet discriminant for the Lax pair [12]. 
Unlike scalar nonlinear Schrodinger equation, the Lax pair for the integrable vector 
nonlinear Schrodinger equations, does not have a complete isospectral theory built 
upon a Floquet discriminant. By restricting uj in the interval (1,2), we are able 
to deal with one unstable mode problem for which only one Melnikov function is 
needed. The Melnikov function will be built out of the Hamiltonian and the invariant 
L 2 norms. 

3. In the 4D invariant subspace of functions independent of x, the dynamics is quite 
different from that of perturbed scalar nonlinear Schrodinger equation [11]. The 
eigenvalues of the saddle are: a positive one of order y/e, a negative one of order 
y/e, and a complex conjugate pair for which the real part is negative, and both real 
and imaginary parts are of order e. The two eigenvalues of order y/e, correspond to a 
fish-shape dynamics, and the other two eigenvalues correspond to a stable spiraling 
dynamics. 

4. The dynamics in the invariant diagonal subspace given by p = q, is governed by the 
PNLS. The unstable manifold of the saddle resides in the invariant subspace. This 
leads to that the persistent homoclinic orbit for the PNLS is the only one for the 
PVNLS, asymptotic to the saddle. 
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2 Dynamics in the 4D Invariant Subspace of Func- 
tions Independent of x 

The 4D subspace 

n = {(p,g) I d xP = d x q = 0} , (2.1) 
is invariant under the PVNLS flow (1.3,1.4). Dynamics in II is governed by 

ipt + \[(\P\ 2 + \q\ 2 ) - uj 2 ] P = te[-ap - (3} , (2.2) 

iQt + \\{\p? + \<l\ 2 ) ~ u 2 }q = ie[-aq - (3} , (2.3) 

2.1 Unperturbed Dynamics 

When 5 = 0, the unperturbed dynamics is governed by 

m + h\p\ 2 + \q\ 2 )-u 2 ] P = o , (2.4) 



2' 

iq t + l[(\p\ 2 + \q\ 2 )-u; 2 } q = 0, (2.5) 



Introducing the polar coordinates 

p = Jhe^ , q = JT 2 e ie > , (2.6) 
one gets the expression for the general solution, 

l\ = constant , I 2 = constant , 

01 = \{h +h ~ w 2 )t + 7i , 

2 = i(/i + / 2 -^ 2 )t + 72 • 
Equilibria consist of the origin p = q = and the submanifold 

s = {(h,i 2 ,e 1 ,e 2 ) \i 1 + i 2 = uj 2 }. (2.7) 

For an illustration of S, see Figure 2.1. 
2.2 Perturbed Dynamics 

In terms of the polar coordinates (2.6), Equations (2.2,2.3) can be rewritten as 



h = e[-2ah -2/tyiiCOS0i] , (2.8) 
0i = \[h + h-u 2 }+e^, (2.9) 



I 2 = e[-2aI 2 -2^I 2 cos6 2 ] , (2.10) 
02 = \[h + h - co 2 ] + e(3^ . (2.11) 



Figure 2.1: An illustration of the equilibrium manifold S. 



Lemma 2.1 All the fixed points of these equations (2.8-2.11) satisfy the conditions 

Ii — I 2 , 6i = 6 2 . 

Proof: Setting the right hand side of (2.8-2.11), one can deduce that for fixed points, 

tan 6i = tan 9 2 , cos 9\ cos 9 2 > , sin 9\ sin 9 2 > . 

Thus 9i — 9 2 , which implies I± — I 2 . □ 

As in the case of perturbed scalar nonlinear Schrodinger equation [11], there are 
fixed points, 

I. The focus £ in the neighborhood of origin, 

a I — 

, cos#i = — — yli , /3 sin 6*i > . 



three 



The eigenvalues are 



. / e(3sin9i\ 2 /— . 

A*i,2 = sa±i* I -j= J -2e(Jyl 1 sm0 1 , 

. (3 sin Q\ 

/1 3A = SOL ± IE j=- , 



where I\ and 9\ are given in (2.12). 
2. The focus P £ in the neighborhood of I± — I 2 — |a> 2 , 

a i — 

= -W 7 i > /5sin0!<O 



' 2 1 eu- 1 J2fP-a 2 u 2 + 



The eigenvalues are 



cos 6*i = 



1 



. /3 sin ^! 
= — ea ± ?e — , 
v A 



„ „ rr . „ ( sin 9i 
■2(3yh sin^x + e I 



where l\ and 9\ are given in (2.15) 



(2.12) 

(2.13) 
(2.14) 



(2.15) 

(2.16) 
(2.17) 



3. The saddle Q £ in the neighborhood of I\ — I 2 — \u 2 , 



h = \uo 2 - euj~ x ^2{3 2 - a 2 uj 2 + ■■■ , cos 
The eigenvalues are 

/i 1)2 = — ea ± \fe 

/i 3)4 = — ea ± is 



h , (3sin9 1 > . (2.18) 



\ 



2(3 y/h sin X -e 
/3 sin 



'/3 sin 6*i 



(2.19) 
(2.20) 



where ii and 9\ are given in (2.18). 

The pair of homoclinic orbits to be located, are asymptotic to the saddle Q £ . Inside 
the subspace IT, attention will be focused upon the neighborhood of I\ = I 2 = \oo 2 . Then 
naturally one wants to introduce coordinates (Ji, J 2 ) in this neighborhood, 



and Equations (2.8-2.11) can be rewritten as 

ii = 



9 1 = 



-auj - 2a Ji - 2(3\l -uj 2 + Ji cos 9 X 



1 / t t x ^ sin #1 

2« + ^7pr^ 



J 2 
#2 



= e 



-au 2 — 2aJ 2 — 2(3\ \uj 2 + J 2 cos 2 



1 / -r -r x ^ sin 6*2 



y/& 2 + h 



(2.21) 

(2.22) 
(2.23) 

(2.24) 
(2.25) 



Rescaling the variables as follows, one can reveal the dynamics at different scales better. 
Let 



Ji = V^ji , -h = Vej2 , r = y/et , 
and Equations (2.22-2.25) can be rewritten as 



f 2 



9' 



where ' = 4-. 

dr 



-au 2 - V2f3u cos 9 l - y/e2aj l - 2(3 \ \\ -uj 2 + v^ji 
^(ji+J 2 ) + V^P^j 



UJ 

V2, 



(2.26) 

cos#i , (2.27) 
(2.28) 



au 2 - V2[3ujcos9 2 - ^2aj 2 - 2(5 (J^uj 2 + y/ij 2 - -^=j cos9 2 , (2.29) 



1 , . . x i- n sin 9 2 

1 yi u + v^J2 



(2.30) 



2.2.1 Leading Order Dynamics 

Only keeping the leading order terms on the right hand side of (2.27-2.30), one gets 



—au 2 — V2(3uj cos 6>! , 



0[ = -U1+J2), 
i 2 = —auj 2 — V2[3uj cos 62 , 



J-2 



9' 



2^1 +h) ■ 



This system is an integrable Hamiltonian system with the Hamiltonian 

H = + j 2 ) 2 + auo 2 {9 1 + 2 ) + V2(3u{smO l + sin# 2 ) 

and another functionally independent constant of motion, 

K = 9 1 -9 2 . 
The Hamiltonian form of (2.31-2.34) is 

OH 



(2.31) 
(2.32) 
(2.33) 
(2.34) 



(2.35) 



(2.36) 



,__&h ,_m 

h ~ 89, ' Ul ~ dj, ' 32 



09 2 



032 



To this order, four lines in the equilibrium manifold £ (2.7) persist. These four lines of 
fixed points of (2.31-2.34) are 



Ji 



'h , 9\ — ± arccos 



9-2 



± arccos 



au 



(2.37) 



Two lines of fixed points for which 9\ 7^ 9 2 will disappear under the full perturbed flow 
(2.8-2.11). One point on each of the other two lines will persist and becomes P e (2.15) or 
Q £ (2.18). The eigenvalues of the four lines of fixed points are 



/xi, 2 = ±y —^/3uj(sm9 1 + sm9 2 ) , ^3,4 = 0. 



(2.38) 



Denote by l u the line in (2.37), for which 9\ = 9 2 and /^sin^x > 0. For l u , n 12 in (2.38) 
are real. Denote by l c the line in (2.37), for which 9\ = 9 2 and /Ssin^i < 0. For l c , /j,i j2 
in (2.38) are purely imaginary. Denote by /q and Iq the rest two lines in (2.37), for which 
9i = —9 2 . For Zq and Zq, // 1j2 in (2.38) are zero. Denote by Q the point on l u such that 
ji = —j 2 = 0. Qq has the coordinates 



3i 



-32 = 0, 9 1 = 9 2 



it — arctan 
— arctan 



y / 2(3 2 -a 



aw 



, if p > , 
, if (3 < . 



(2.39) 
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Figure 2.2: The stable and unstable manifolds of I. 



Figure 2.3: The center manifold of l u . 
Denote by Po the point on l c such that j\ = —ji = 0. Po has the coordinates 

'y/2/3 2 -a 2 i. 



h = -32 = 0, 9 1 = e 



2 



7r + arctan 



arctan 



, if < . 



, if p > , 

(2.40) 



Next we study the dynamics of (2.31-2.34). Denote by Ha the 3D subspace of IT (2.1), 

n A = {(ji,j 2 ,#iA) I if = 0i -0 2 = A} . (2.41) 

n A is an invariant subspace under the flow (2.31-2.34). Restricted to IT A , H (2.35) takes 
the form similar to that in [11], 

H = ij 2 + 2au 2 9 + 2V2f3tu cos (A/2) sin0 , 

where j = j x + j 2 and Q = Q x - A/2. Denote by W^(l u ), Wfi(l v ), and Wg(Z„) the 
2D stable, unstable, and center manifolds of l u respectively. W^(l u ) = Wjj(l u ) has the 
topology shown in Figure 2.2. W^(l u ) has the topology shown in Figure 2.3. Similar to 
W^(l u ), there is also an invariant 2D submanifold W^{l c ) around l c as shown in Figure 
2.4. W^(l u ) = Wfi(l u ) is given by the singular level set of H in n A =o- In n A (A ^ 0), 
similar picture holds as shown in Figure 2.5, but with drifting along the asymptotic line. 
Denote by W^(l u ) and W^{l u ) the 3D center-stable and center-unstable manifolds of l u . 
Wn(l u ) = Wyi{Iv) has the topology shown in Figure 2.6. 

2.2.2 Full Dynamics 

Since W^{l u ) is normally hyperbolic under the leading order flow (2.31-2.34), it persists 
under the full perturbed flow (2.27-2.30). Denote this persistent 2D locally invariant 
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Figure 2.4: The center manifold of l c . 



Figure 2.5: Stable and unstable manifolds in Ha (A ^ 0). 



Figure 2.6: The center-stable and center-unstable manifolds of I 



Figure 2.7: Persistent center-stable and center-unstable manifolds. 
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center manifold by (l u ). Denote the persistent 3D locally invariant center-stable 

and center-unstable manifolds by W^ s ' £ \l u ) and W^ £ \l u ). W^ s ' £ \l u ) and W^' e \l u ) 
have the topology shown in Figure 2.7. Notice that W^ s,£ \l u ) forms the stable manifold 
of the saddle Q e (2.18), which we also denote by Wf[(Q £ ). The unstable manifold of Q £ , 
denoted by W^(Q £ ), is one dimensional. 
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3 Global Integrable Theory 

Consider the integrable vector nonlinear Schrodinger equations (e = in (1.3,1.4)), 

ipt + Pxx + ^[(\p\ 2 + \q\ 2 ) ~uJ 2 ]p = , 



iqt + q xx + ^ [(b| 2 + \q\ 2 ) ~ " 2 ]<1 = • 



The Lax pair is given as [13], 



where 



U = \A + A-y , V = \ 2 A + \A ± + A 2 , 





( 


-¥ 





^ 


(° 


\v 




\ 


A = 







H 





> Ai =\ -\P 












V 








¥ ) 










/ 



A, = - 



/ -(\p\ 2 + \q\ 2 ) + Alu 2 -2p x -2q x 
-2p x \p\ 2 + 2u 2 qp 

V -2q x qp \q\ 2 + 2u 2 



3.1 Linearized Equations 

We start with the solution of (2.4,2.5), 

p = ae iSl , q = be i52 , 

where 

Si = -(a 2 + b 2 - uj 2 )t + 7l , 5 2 = -(a 2 + b 2 - u 2 )t + 72 . 
Linearizing (3.1,3.2) at this solution, by setting 

p = e i5 i[a+p] , q = e^[b + q] , 

one gets 



ipt + Pxx + 2 a [ a (p + p) + Kq + <?)] = » 



Setting 



«<Zt + + t^Kp + P) + % + ?)] = • 



p = f +e ikx+nt + j_ e -ikx+cit ^ - = g+e ikx+nt + _ e - lfca; +™ ? (k eZ ) 



(3.1) 
(3.2) 

(3.3) 



(3.4) 



one gets the dispersion relations 



fii, 2 = ±£Va 2 + 6 2 - k 2 , fi 3 ,4 = ±^ 2 



(3.5) 
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Thus, if we choose 1 < \Ja 2 + b 2 < 2, then there is only one unstable mode cosrr. This is 
the reason that we restricts u as in (1.1,1.2). When k — 1, 

Q lj2 = ±<r , a = Va 2 + b 2 -l . (3.6) 

The corresponding eigenfunctions are 

p = ae ±llp e ±at cosx , q = be ±llp e ±ut cosx , <p = arctanc . 



3.2 Homoclinic Orbits and Figure Eight Structures 

The Backlund-Darboux transformation given in [16] [6] will be utilized to construct ho- 
moclinic orbits asymptotic to the periodic orbits (3.4) up to phase shifts. 

Theorem 3.1 ([16] [6]) Let p and q be a solution of the vector nonlinear Schrddinger 
equations (3.1,3.2), and let tp = (ipi, tp2, i>3) T be an eigenfunction of the Lax pair (3.3) at 
(\,p,q). If one defines 

P = ^A-A W+ g +fer ,3.7) 

* = « + *<*- A W + KiW - (3 ' 8) 

then p and q also solve the vector nonlinear Schrddinger equations (3.1,3.2). 

We solve the Lax pair (3.3) at (A,p, q) with p and q given by the periodic orbits (3.4), we 
get 

( i>f \ ( (i\Ti)e i(5l+S2) \ 



4 



± 



ae 



\if>f J V be^ J 



e ±ait+i<T2t e ±iK ± x 



where 



« ± = -^A± lVa 2 + b 2 + X 2 , Va 2 + b 2 + A 2 = 1 , (3.10) 
6 2 

*i = y, a 2 = -l[7(a 2 + 6 2 ) + 2]. (3.11) 

Remark 3.1 From (3.11), we see that in order to have temporal growth, the imaginary 
part of A can not be zero. Also in order to have a nontrivial Backlund-Darboux transfor- 
mation (Theorem 3.1), the imaginary part of\ can not be zero. From (3.10), ifa 2 +b 2 > 1, 
then A is purely imaginary. If we also require a 2 + b 2 < 4, i.e., in the one unstable mode 
regime (3.5), then Va 2 + b 2 + A 2 > 2 will lead to real A and no temporal growth. 

Remark 3.2 Notice that the imaginary part of is not zero, thus ■0 ± are not peri- 
odic or antiperiodic functions in x. In fact, they grow or decay in x. This fact shows 
that complex double points are not necessary in constructing homoclinic orbits through 
Backlund-Darboux transformations (see also [9]). 
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Let 



ip = c ip + c ip , 



where c + and c are two arbitrary complex constants. We introduce the notations 

c +/ c - = e ^ , A = -ia , a= Va 2 + b 2 - 1 , 
a + j = Va 2 + &V ( ^ o) , = arctan a , a x = a/2 , 



a p Id 

r = 7r^+o? y — n x + 



2 ' 



(3.12) 

(3.13) 
(3.14) 

(3.15) 



where p and $ are called the Backlund parameters. If one chooses A = ia, one will get 
the same result. We can rewrite ip in the following form, 



^3 



7T 7T 

cosh r cos(y h i? ) + * sm h t sin(y h i? ) 

2 2 



2Vc+c- Va 2 + 62 e i(«i+«2) e i^* e -^x ? 



cosh r cosy + isinhrsiny 2v / c" t "c ae""' z e""' z "e 6 " 
cosh r cos y + i sinh r sin y 2V c+c~ be tSl e l(T2t e~^ ax , 

where the factor 2\ / c + c~e la2t e~^ ax will be canceled away. By (3.7,3.8), one gets 

p = ae iSl h, q = be iS2 h, 

where 



(3.16) 



h 



cos(2^ ) + i sm(2# ) tanh(2r) — sin-# sech(2r) cos(2y + # — 



7I\ 



7I\ 



1 + sin-# sech(2r) cos(2y + $ ) 

2 



(3.17) 



where 5\ and 5 2 are given in (3.4), other notations are given in (3.13-3.15). For the 
moment, the even-in- a; condition has not been inforced. For fixed a and b, the phases 
7x and 72 (3.4), and the Backlund parameters p and d parametrize a four dimensional 
submanifold with a figure eight structure. If a = b and 'ji = 7 2 , then the homoclinic orbit 
reduces to that for the scalar nonlinear Schrodinger equation [10]. 
As t — > ±oo, 



h 



o ±i2$ 



(3.18) 



The even-in- x condition can be achieved by restricting the Backlund parameter d to 
special values. That is, if one requires that 



TX 



then h is even in x, 



h 



# + A) - - = , 71 , 



cos(2-# ) + ?sin(2^ ) tanh(2r) =F sin^ sech(2r) cos a; 



(3.19) 



1 ± sin-^o sech(2r) cosx 



(3.20) 
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where the upper sign corresponds to d + $0 — f = 0. For fixed a and b, the phases 71 
and 72 (3.4), and the Backlund parameter p parametrize a three dimensional submanifold 
with a figure eight structure. If one further chooses 71 = 72 (3.4), then (3.4) represents a 
periodic orbit. In such case, the phase 71 = 72 and the Backlund parameter p parametrize 
a two dimensional submanifold with a figure eight structure. 



3.3 A Melnikov Vector 



The Hamiltonian for the integrable vector nonlinear Schrodinger equations (3.1,3.2) is 
given as, 

H= {\p x \ 2 + \q x \ 2 - i[(|p| 2 + |g| 2 ) 2 - 2c 2 (b| 2 + |g| 2 )]bi* . (3.21) 
Jo 4 



(3.1,3.2) can be rewritten in the Hamiltonian form, 



5H 



5H 



5H 



5H 



iPt 



-IF" » % Vt = --r- , iqt 
op op 



Sq ' ^ Sq 



The two L 2 norms 



t>27T l>27T 

E 1 = \p\ 2 dx , E 2 — / |g| 2 ofe , 
io Jo 



(3.22) 



are also invariant functionals. We will build the Melnikov vector by the gradient of the 
invariant functional G obtained through a combination out of H, E±, and E2, 



G = H + 



47T 



+ e 2 ) - ^ 



(E 1 + E 2 ) 



f*{\p*\ 2 + \Q*\ 2 -\(\P\ 2 + l^l 2 ) 2 }^ + ^(^1 + E i? ■ (3-23) 



16 



Figure 4.1: The graph of a — ol{uj). 

4 Persistent Homoclinic Orbits 

4.1 The Invariant Diagonal Subspace 

The diagonal subspace 

E = {(p,q) \p = q} 

is invariant under the PVNLS flow (1.3)-(1.4). The dynamics in E is governed by the 
PNLS 

ift + + ^PM 2 ~ = ie ferx - - /3] , (4.1) 

studied in [10]. 

4.2 The Main Theorem 

Theorem 4.1 (Main Theorem) There exists a e > 0, snc/i that for any e G (0, £o)> 
there exists a codimension 1 surface in the space of (a, [3, uj) G R + x R + x R + where 
uj G (1,2)/ S, S is a finite subset, and acu < \f2j3. For any (a,(3,uj) on the codimen- 
sion 1 surface, the singularly perturbed vector nonlinear Schrodinger equations (1.3)-(1.4) 
possesses a homoclinic orbit asymptotic to the saddle Q € (2.18). This orbit is also the 
homoclinic orbit for the singularly perturbed scalar nonlinear Schrodinger equation (4-1) 
studied in [10], and is the only one asymptotic to the saddle Q e (2.18) for the singularly 
perturbed vector nonlinear Schrodinger equations (1.3)-(1.4). The codimension 1 surface 
has the approximate representation a = a(cu) given in Figure 4-1, obtained in [10]. 

Proof. Existence of a homoclinic orbit has been proved in [10]. Under the PVNLS 
flow (1.3)-(1.4), the saddle Q t (2.18) has a two dimensional unstable manifold W U (Q £ ) 
which is included in the invariant diagonal subspace E, 

W U (Q £ ) C E . 

Thus any homoclinic orbit asymptotic to Q e has to be inside E. □ 

Remark 4.1 Following the same arguments as in [10], Unstable Fiber Theorem and 
Center-Stable Manifold Theorem can be easily established for the PVNLS (1.3)-(1.4) ■ 
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4.3 An Alternative Melnikov Function 



In [10], the Melnikov function was built with the Floquet discriminant. Here we give an 
alternative Melnikov function built with the invariant functional G: 



M 



r -dt = r i 

J-oo dt J-oo JO 



oo r2w 
oo J 



5G 8G_ 8G 8G_ 

-r-Pt + -jzPt + -r-qt + -jrqt 
op op oq oq 



dxdt , 



where the integral is evaluated along the homoclinic orbit (3.16,3.20) at a = b = -^u and 
5\ = 71 = 7 2 = 5 2 , ± in (3.20) leads to the same result. 



M = e 



where 



Mi = -uj 
2 



M 2 = tu : 



1 r+oo r 2iT 

L 2 



oo JO 

oo r2w ' 



Mi - aM 2 - p2V2u cos 71 M 3 



u 2 (l - \h\ 2 )(hh xx + hh xx ) - A\h x 



(4.2) 



: r°° r 

J-00 Jo 



M 3 



uj 2 \h\ 2 {l-\h\ 2 )-{hh xx + hh xx ) 
00 r2-K r\ 



dxdr 
dxdr , 



/ + OO rZTT \ , s 

-00 /o [r 2 h {r \l-\h\ 2 )-h[ 



r) 

XX 



dxdr , 



where is the real part of h (3.20). 

In [10], we have showed that M = together with 

H(0, 0, 71 + 2^o) - H(0, 0, 71 - 2^ ) = , 



(4.3) 



give the function a = a(u>) which predicts the leading order location of the homoclinic 
orbit in the external parameter space. Here ±2-# are the asymptotic phases (3.18), and 
the Hamiltonian 7i is given in (2.35), 

n{ji,j 2 ,di) = ^(Ji+h) 2 + 2c^ 2 #i + 2v / 2/3wsiii0i . 



(4.3) leads to 

COS71 = 

Combining (4.4) with M = 0, one gets 

a = a (u) = Mi M 



(3 sin 2^ ' 
sin 2-^0 



M s 



(4.4) 



(4.5) 



which has the same graph as in Figure 4.1. Thus the Melnikov function above is equivalent 
to that in [10]. 
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5 Conclusion and Discussion 



The instability of the singularly perturbed vector nonlinear Schrodinger equations (1.3)- 
(1.4) happens along the diagonal direction {p = q), which leads to the fact that the 
homoclinic orbit in the invariant diagonal subspace {p = q) is the only one for (1.3)-(1.4). 
Existence of a homoclinic orbit in the invariant diagonal subspace (p = q) governed by 
the singularly perturbed scalar nonlinear Schrodinger equation (4.1) has been established 
in [10]. 

If one releases the even constraint p(t, —x) = p(t, x) and q(t, —x) = q(t, x), the homo- 
clinic orbit is given by (3.17) in which the spatial variable x appears in the combination: 

x + + A) - | , 

where d is the Backlund phase-parameter. (Even constraint is realized by setting d + 
t? — | = or 7r.) i? characterizes the x-location of the top of the hump h(t,x) for any 
t. It is clear that the Melnikov integrals evaluated along h(t, x) are independent of d. In 
fact, since the PVNLS (1.3,1.4) is translationally invariant in x, the persistent homoclinic 
orbit is not isolated, if it exists. In the current study, the difficulty associated with 
noneveness, that we are facing, is that we need two Melnikov integrals due to the extra 
phase-parameter Using the real and imaginary parts of the constant of motion F\ (see 
(3.22) in [10]), we can build two Melnikov integrals. Unfortunately, we found numerically 
that the second Melnikov integral built with the imaginary part of i*\ is identically zero. 
We do not know whether or not one can build two non-identically-zero Melnikov integrals. 
If not, second or higher order Melnikov integral might be needed. 

The matter is much more important than just a technical difficulty. In this noneven 
case, Ablowitz et al. [2] [1] found a novel scenario in which chaos occurs, that is, chaotic 
flipping between left and right travelling waves. This scenario is generic with respect to 
perturbations [2] [1]. When even constraint is not imposed, this new scenario of chaotic 
dynamics is observed in experiments [1]. With the difficulty mentioned in the last para- 
graph and the new results of Ablowitz et al. [2] [1], understanding the chaotic dynamics 
in the noneven case becomes a much more challenging and significant problem. 



Acknowledgement: The author is greatly indebted to Professor Mark Ablowitz and 
Professor Jianke Yang for fruitful discussions which led to the current study. 

Figure Captions: 

Figure 2.1: An illustration of the equilibrium manifold £. 

Figure 2.2: The stable and unstable manifolds of l u . 

Figure 2.3: The center manifold of l u . 

Figure 2.4: The center manifold of l c . 

Figure 2.5: Stable and unstable manifolds in Ua (A^O). 

Figure 2.6: The center-stable and center-unstable manifolds of l u . 

Figure 2.7: Persistent center-stable and center-unstable manifolds. 

Figure 4.1: The graph of a — a (u). 



19 



References 

[1] M. J. Ablowitz, J. Hammack, D. Henderson, and C. Schober. Modulated Periodic 
Stokes Waves in Deep Water. Physical Review Letters, 84, no. 5:887-890, 2000. 

[2] M. J. Ablowitz, B. Herbst, and C. Schober. Computational Chaos in the Nonlinear 
Schrodinger Equation Without Homoclinic Crossings. Physica A, 228:212-235, 1996. 

[3] M. J. Ablowitz, Y. Ohta, and A. D. Trubatch. On Discretizations of the Vector 
Nonlinear Schrodinger Equation. Phys. Lett. A, 253, no. 5-6:287-304, 1999. 

[4] M. J. Ablowitz, Y. Ohta, and A. D. Trubatch. On Integrability and Chaos in Discrete 
Systems. Chaos Solitons Fractals, 11, 1-3:159-169, 2000. 

[5] M. G. Forest, D. W. McLaughlin, D. J. Muraki, and O. C. Wright. Nonfocusing 
Instabilities in Coupled, Integrable Nonlinear Schrodinger pdes. J. Nonlinear Sci., 
10, no.3:291-331, 2000. 

[6] M. G. Forest, S. P. Sheu, and O. C. Wright. On the Construction of Orbits Homoclinic 
to Plane Waves in Integrable Coupled Nonlinear Schrodinger System. Phys. Lett. A, 
266, no. 1:24-33, 2000. 

[7] A. Hasegawa and Y. Kodama. Solitons in Optical Communications. Academic Press, 
San Diego, 1995. 

[8] M. N. Islam. Ultrafast Fiber Switching Devices and Systems. Cambridge University 
Press, New York, 1992. 

[9] Y. Li. Backlund-Darboux Transformations and Melnikov Analysis for Davey- 
Stewartson II Equations. Journal of Nonlinear Sciences, 10, no. 1:103-131, 2000. 

[10] Y. Li. Persistent Homoclinic Orbits for Nonlinear Schodinger Equation Under Sin- 
gular Perturbation. Submitted, 2001. 

[11] Y. Li, D. McLaughlin, J. Shatah, and S. Wiggins. Persistent Homoclinic Orbits for a 
Perturbed Nonlinear Schrodinger equation. Comm. Pure Appl. Math., XLIX:1175- 
1255, 1996. 

[12] Y. Li and D. W. McLaughlin. Morse and Melnikov Functions for NLS Pde's. Com- 
mun. Math. Phys., 162:175-214, 1994. 

[13] S. V. Manakov. On the Theory of Two-Dimensional Stationary Self-Focusing of 
Electromagnetic Waves. Soviet Physics JETP, 38, no.2:248-253, 1974. 

[14] C. R. Menyuk. Nonlinear Pulse Propagation in Birefringent Optical Fibers. IEEE 
J. Quantum Electron, 23, no.2:174-176, 1987. 

[15] C. R. Menyuk. Pulse Propagation in an Elliptically Birefringent Kerr Medium. IEEE 
J. Quantum Electron, 25, no.l2:2674-2682, 1989. 



20 



[16] O. C. Wright and M. G. Forest. On the Backlund-Gauge Transformation and Homo- 
clinic Orbits of a Coupled Nonlinear Schrodinger System. Phys. D, 141, no. 1-2: 104- 
116, 2000. 

[17] J. Yang and Y. Tan. Fractal Dependence of Vector-Soliton Collisions in Birefringent 
Fibers. Phys. Lett. A, 280:129-138, 2001. 



21 



